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O , We present a multilateral theoretical study of bound polarons in oxide compounds MgO and 

a-Al2 03 (corundum). A continuum theory at arbitrary electron-phonon coupling is used for cal- 
K^ culation of the energies of thermal dissociation, photoionization (optically induced release of an 

Q^ ■ electron (hole) from the ground self-consistent state), as well as optical absorption to the non- 

t^ . relaxed excited states. Unlike the case of free strong-coupling polarons, where the ratio k, of the 

f>^ \ photoionization energy to the thermal dissociation energy was shown to be always equal to 3, 

^D ' here this ratio depends on the Frohlich coupling constant a and the screened Coulomb interaction 

strength (3. Reasonable variation of these two parameters has demonstrated that the magnitude 
of K remains usually in the narrow interval from 1 to 2.5. This is in agreement with atomistic 
calculations and experimental data for hole 0~ polarons bound to the cation vacancy in MgO. 
The thermal dissociation energy for the ground self-consistent state and the energy of the optically 



fi 

'^ \ induced charge transfer process (hops of a hole between 0^~ ions) have been calculated using the 

Q ' quantum-chemical method INDO. Results obtained within the two approaches for hole 0~ po- 

O ■ larons bound by the cation vacancies (V~) in MgO and by the Mg^"*" impurity (V Mg) in corundum 
are compared to experimental data and to each other. We discuss a surprising closeness of the 

KA , results obtained on the basis of independent models and their agreement with experiment. 

c5 ■ PACS numbers: 73.61.Ng,71.38.-hi,71.15.-m 
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Properties of polarons and bipolarons in oxide materials continue to attract considerable 
attention in solid state and materials physics due to their possible relevance to the high- 
Tc superconductivity |^. Theory for large-radius [Q] and small-radius polarons is well 
developed. At moderate temperatures polarons are usually bound to defects (vacancies or 
impurities) unavoidably present in any material. Examples are a hole trapped by the cation 
vacancy in MgO (called V~ center) and a hole trapped by the Mg^+ impurity substituting 
for a regular Al^+ atom in AI2O3 (called Ymq center) ^. 

Fundamental characteristics of free polarons are energies of their photoionization (-Eopt) 
and thermal dissociation {El^. From the well-known theorem 1:2:3:4 derived by Pekar ||^ for 
relations between the average kinetic energy of the electron, the polarization energy of the 
crystal -Epoiariz' ^^^ modulus of the electron eigenenergy El, and the average potential energy 
of the interaction between the electron and the polarized crystal in the strong-coupling 
polarons, it follows, in particular, that the ratio between the photoionization energy E^^^ = 
\El\ and the thermal dissociation energy Ei^ = \El\ — -Epoiariz ^^ ^ fr^e polaron k-^ = 
EJ^^^/eI^ = \EI\/{\EI\ — -Epoiariz) equals exactly 3. A generalization of this ground-state 
theorem, which relates the interaction energy to the kinetic energy for free polarons, was 
derived by Lemmens, De Sitter and Devreese for all values of the electron-phonon coupling 
constant 0]. Optical transitions of large polarons with taking relaxed excited states into 
account were examined in Refs. [0]. 

Calculation of the ratio k = Eopt/Eth of the photoionization and thermal dissociation 
energies allows one to examine the adequacy of different theoretical approaches to bound 
polarons. In the present Letter we analyze this problem in two ways: using an extended 
continuum theory of large radius bound polarons, and performing atomistic calculations for 
the particular hole polarons in the above-mentioned ionic solids — MgO and corundum 
(a;-Al203). These two types of calculations allow us to draw conclusions how the continuum 
approach can describe models, where the lattice polarization is of Frohlich type, while the 
state of charge carriers is of small radius (such a model has been suggested e. g. in Ref. |@] 
for a small Frohlich polaron). In the framework of this approach, unknown values of the 
macroscopic parameters should be obtained from comparison with experimental data. At 
the same time, we examine the efficiency of the microscopic quantum-mechanical approaches 
for solving the aforementioned problems of bound states. 

A variational theory of bound polaron states which embraces the overall interval of values 
of the Frohlich coupling constant a and of the strength (3 of the Coulomb interaction between 
a band charge carrier and a charged impurity was developed in Ref. Q. The Hamiltonian of 
an electron (hole) in the field of an impurity interacting with the polar longitudinal optical 
(LO) phonons is 

2 2 

i^ = TT- - — + E ^^Lo4«k + E(^A:ake^'^" + h.c), (1) 

where r is the position coordinate operator of the electron (hole) with band mass m^, p 
is its canonically conjugate momentum operator; a\ and ak are the creation (and annihi- 
lation) operators for LO phonons of wave vector k and frequency cUlo- The Vk are Fourier 
components of the electron(hole)-phonon interaction 
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(2) 
with the Frohhch couphng constant 

a = -^^ (^ - 1) ,/^, (3) 

where Eq and £oo are the static and high-frequency dielectric constants of the crystal, respec- 



tively. Choosing the Feynman units, when the polaron radius Jh/{2mi,Ui^o) and the phonon 
energy ^cUlo serve as units of length and energy, respectively, we represent the second term 
in the right-hand side of Eq. (|l]) in the form —f3/r, where the dimensionless Coulomb 
interaction strength relates to the Frohlich coupling constant (^ as follows: 

P = 2a/ieo/e^ - 1). (4) 

The approach, developed by Devreese et al. 0, consists in two stages. First, the unitary 
transformation is performed 



U = exp 



Yl (^fcflkPk - h.c.) 



Pk = {(Pn\e'n(pn), (5) 



where |0„) is the electron (hole) wave function, which will be found variationally. This 
transformation of the Hamiltonian (|I|) results in 

H' = U-'HU = Ho + J2 4«k + Hi^t, (6) 

k 

i^o = P^ - - + E W\P^\' - E iVkl'ipt^''-'' + c.c), (7) 

"^ k k 

Hint = E(^fc«k(e*''" - Pk) + h.c). (8) 

k 

Second, the following trial function is proposed in Ref. §], which is proved to be a fair 
approximation for all values of a and (3: 

1^) = c|O)|0„) + E n*^k(e~'''-'' - Pimm, (9) 

k 

where c is a normalization constant and |0) is the phonon ground state. The variational 
principle with the trial function (^ results in 

9^ 1 - W? (10) 



where 



c A(k) + D2(k) + (l-|pkP)x/2' 



X=-El^fcl'(l-|PkP)(^k + ^k), (11) 
^ k 

D^{k) = (0„|(e^^-^ - Pk)i^o(e-^'^-'- - plMn), (12) 

D^ik) = (1 - |pkn(l - {K\Ho\(pn)). (13) 



The variational energy for the bound polaron is found to be 

E={<j)n\Ho\<Pn) + ^. (14) 

This energy is to be minimized with respect to the parameters entering the trial function 
|0„). It is this minimal energy which relates to the thermal dissociation of a polaron bound 
to an impurity: Eth = |-Emin|- The electron (hole) eigenenergy is Ec = -Emin — -Epoiariz, where 
-S'poiariz is the elastic energy of the polarization (i.e., the energy of the ionic polarization 
due to a lattice energy gain after self-consistent relaxation). The minimal energy, which is 
absorbed during the photoionization of a bound polaron (the optically induced release of an 
electron (hole) from the ground state), equals the modulus of the electron (hole) eigenenergy 
\Ec\: 

Eopt = \Ec\ = Eth + -Epolariz- (15) 

Hence, the ratio of the photoionization energy and the thermal dissociation energy is 

-'-'opt -C'th ~r -C/polariz /■\n\ 

^ = ^^ = ^ • (16) 

In order to find the ratio k, it is necessary to calculate the eigenenergy and the energy 
of the polarization field. Under the unitary transformation of Eq. (^), the Hamiltonian of 
the phonon subsystem turns to 

Kh = U'' E 4«kf/ = E 4«k + E l^^riPkr - E(^^Pkak + h.c). (17) 

k k k k 

Averaging the transformed Hamiltonian of the phonon subsystem on the trial function (^ 
gives the energy of the polarization field in a polaron bound to an impurity: 

^polariz = imUl^) = Y^ + E I^^HPkP (18) 

with 

IT/ \2(^ I „ |2\3 

^ [D,{k) + DM - {1 - W\^)x/2f ^ ^ 

The variational parameters which deliver a minimum to the polaron energy, should be sub- 
stituted to Eq. (|1^). In particular, for a spherically symmetrical variational function |0„) 
the expressions in the right-hand side of Eq. (0) can be simplified to 

x='-^r dk ^i^Mi ^, (20) 

ElHnPkP = -rrfA:|p,p. (21) 

k ^ ^0 

Using for the ground state the hydrogen-like function, which was suggested in Ref. M, 



1^..) = ^r"' (22) 

with a variational parameter a, we obtain the ratio K.{a,l3) = Eopt,is/Eth,is as a function of 
the parameters a and (3 for a bound large polaron. 

Results of the numerical calculation of the ratio K{a,l3) for a polaron bound to an 
impurity at < a < 10 and < (3 < 10 are plotted in Fig. 1. As follows from this figure, 
for all a 7^ the ratio K{a,P) > 1. This inequality arises due to the following fact. During 
the photoionization process, the electron (hole) is released from the polarization potential 
well which is frozen, being adjusted to the electron (hole) wave function in the polaron 
ground state, while in the case of a thermal dissociation of the bound polaron there occurs a 
disintegration of the polarization state and, consequently, the coupling of the electron (hole) 
to the polarization well weakens continuously during this process. 

As seen from Fig. 1, a decrease of the Coulomb attraction parameter /3 at a fixed value 
of a leads to a rise in the ratio K,{a,l3), except for the region of (3 less than or of the 
order of unity at a < 8. In particular, at a = 10 and j3 = 0, K,{a,P) is as large as 2.6, 
already close to but still smaller than 3, obtained in the polaron theory in the limit a — > cxo. 
When increasing a, at a fixed value of (3, a role of the interaction of an electron (hole) with 
polarization enhances, and thus, the value of the ratio K,{a,l3) rises. A weak maximum in 
the region of small values of /3 ~ 1 at a < 8 apparently reflects a transition from a two- 
center state of a polaron moving around the impurity center to a one-center state, where an 
electron (hole) is in a unique potential well due to the common action of the polarization 
field and of the impurity center. 

Further, we consider the excited state of the bound polaron in the non-relaxed polariza- 
tion well of the Is ground state of Eq. ([l^) , which will be referred to as a non- relaxed excited 
state. The factors Vj^p^. at aj^ and V^pk at a^ in the unitary transformation of Eq. (^ have 
the meaning of shifts of the polarization vibrations due to the interaction of the polarization 
field with the electron in the state </?„: 

U-'a^U = ak - Vj:pI, U~^alU = aj, - \4pk. (23) 

Hence, the wave function of the non-relaxed excited state 2p in the frozen polarization well 
correlated to the electron state Is has the form: 

|^2p) = c|O)|02p) + E V^;^k^k|O)|02p), Ak = e-''^-^ - p*^% 

k 

where p]^ = {4>is\e^^'^\(pis)- Like in the case of self-consistent polaron states, gi^ is a vari- 
ational function. After transformations similar to the above-described ones we obtain the 
following variational energy of the non- relaxed excited 2p state: 



As 
2 
where in 



El; = {hp\Ho\hp) + ^, (24) 



6 is a variational parameter and 6 is the angle between r and the 2;-axis. In Eq. 



;^i^ = -- r dk r sin ede^r—- ^^^^^""^ , ^ . . 1 . , (25) 

Di(k) = (02p|Ak//o^kl02p), (26) 

D2(k) = (A*Ak)(l - (02p|^o|02p)). (27) 

The thermal dissociation of the excited state (2p) in the polarization well adjusted to the 
ground state (Is) is (^th)2p = l^2p,mml- 

The elastic energy of the polarization field is 



(-E^polariz)2p = T" T ^ XI I ^^ I \P^\ ' (^8) 



which is similar to that for the Is state. The first term in the right-hand side of this equation 
is by two orders of magnitude smaller than the second one, therefore in both cases 

^polariz-El^fel'lPkl'- (29) 

k 

The electron (hole) eigenenergy, corresponding to the non-relaxed excited state 2p of a bound 
polaron, is 

(-^c)2p — -^2p,min ^ -C'polariz, (30) 

and the minimal energy, which would be absorbed during the photoionization of a bound 

)2p — \\^c)2-p\^ 



polaron from its non-relaxed excited state 2p, (i?opt)2p = l(-^c)2pl5 takes the form 



(-E'opt)2p — (-^th)2p + -^polariz- (31) 

Calculations for large polarons were performed with the following parameters, (i) For 
MgO, ^cjlo = 0.09eV, £00=2. 96 and £o=9.86 |jlO|, "^5 = 2.77mo lll|, where mo is the bare 



electron mass; according to the definition (^) a. = 4.83, whereas the calculation by Eq. (^ 
gives (3 = 4.14. (ii) For a-A^Og, hui^o = 0.07eV|3, £oo=3.1 and eo=9.0 0. Estimate 
of the hole effective mass in corundum from the experimental value Eopt = 2.58 eV gives 
rrih = 4.0mo, and consequently, a = 5.41; /5 = 5.66. 

Atomistic calculations for bound hole polarons were performed using the quantum chem- 
ical method of the Intermediate Neglect of the Differential Overlap (INDO) [|1^] modified 



for ionic solids ||15[. This is based on the Hartee-Fock formalism and allows self-consistent 
calculations of the atomic and electronic structure of pure and defective crystals, that is 
we are able to simulate polarons making no a priori assumptions about their geometrical 
structure. The INDO method was successfully applied for the study of many defects in oxide 



crystals (see a review article |]T6[ and references therein). Quantum cluster is embedded into 



^This value is suggested in Ref. [|T2J as an average of the phonon energies corresponding to three 
modes which provide the strongest contributions to the dielectric dispersion. This average value 
is in good agreement with the more recent data |jl^, where the measured optical phonon energy 
range from 0.052eV to O.lOSeV. 



the electrostatic field of the infinite non-point crystalline lattice. To model a hole polaron, 
we remove one electron from the cluster and allow all atoms to relax to the minimum of the 
total energy which gives us the atomic structure of a polaron. The analysis of a relevant final 
charge density distribution characterizes the polaron's electronic structure. Calculations of 
the difference of total energies for the ground and excited states with a fixed lattice geometry 
permit, in principle, to find the activation energy Eopt of a charge transfer process, when a 
hole after optical stimulation hops between 0^~ ions. In its turn, the thermal dissociation 
energy _Eth is calculated as the energy gain due to atomic displacements induced by a hole. 
Results of the INDO method are illustrated in Fig. 2. 

In MgO, calculations were done for 125- atom clusters, modelling 9 spheres of atoms 
surrounding a cation vacancy in the coordinate origin. The basis set is the same as in 
our previous MgO studies Ijl^. In corundum, 65-atom stoichiometric clusters containing 



13 formula units were used. The INDO parameters and basis set remain also the same as 
in previous a-Al203 calculations |T3]. Both our microscopic calculations and ESR data P| 



demonstrate that a hole is well localized by a single oxygen atom, forming 0~ ion consider- 
ably displaced from its regular lattice site. Practically, it is difficult to achieve a convergence 
in the self-consistent-field calculations for the delocalized (band) states necessary for calcu- 
lations of the photoionization energy. This is also the case in our calculations for MgO and 
corundum. Neither we are aware of the relevant experimental data on optical ionization 
energies of these polarons. 

Results of calculations for the V~ centers in MgO by theories CTAC, CTSC and INDO 
and the available experimental data are collected in Table I. 

The thermal dissociation energy Eth obtained in the continuum theory is close to the 
experimental value; the difference of about 10% between them seems to be within the limits 
of the error of measurement. The INDO method gives the value of E'th which is also close 
to the experimental data. The energy Eopt calculated in the continuum theory and within 
the INDO approach are both in a good agreement with experiment. It is worth noting that 
CTAC better compares with experiment than CTST. 

A detailed discussion of the optical properties of the V~ centers in MgO has been pre- 
sented in terms of small-radius polarons by Schirmer ||19| . Theory predicts the ratio of the 
optical to thermal reorientation energy to be four, and the ratio to polarization energy to 
be two. The latter relation is also supported by our CTSC calculations (Table I). A very 
qualitative discussion of the large-radius vs small-radius hole polarons in corundum has been 
done in Ref. fM . 



As distinct from the case of MgO, calculation of the polaron energy spectra in corundum 
lacks for some material parameters. The band mass rrif, has not been measured; the value 



huji^o = 0.07eV |T^ is not fully reliable; the thermal dissociation energy for the ground state 
of the Vmq center is not known. The value of 0.7 eV from Table II is, in fact, the activation 



energy of a free small polaron |^, which can be only conventionally considered as E^y^. It 
is worth noting, that the shell-model calculation underestimates the energy Et^, resulting 
with 0.56 eV ||22[. The energy i^th calculated by using the continuum theory is close to 
the respective result of the INDO approach; the both of them differ from the experimental 
value. 

In conclusion, we have shown that the continuum arbitrary-coupling theory (CTAC) of 
large bound polarons qualitatively and quantitatively describes the energy spectrum and 



transition energies for the impurity centers in oxide compounds MgO and Q;-Al203. The 
thermal dissociation energies obtained in the framework of CTAC are rather close to the 
results of INDO. The CTAC has the advantage that it allows to describe the impurity energy 
spectrum in detail. The surprising fact, that two essentially different approaches have led 
to close results, seems to be indicative of a dualistic nature of the bound polarons, which 
combine a macroscopic scale of polarization changes with microscopic (atom-like) charge 
carrier states. The efficiency of the proposed in our letter macroscopic approach to the 
bound polarons in the aforementioned oxides is probably due to the fact, that the specific 
features of the small-radius (atom-like) states of charge carriers do not appreciably influence 
the distribution of the polarization field that is induced by themselves. 
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FIG. 1. Ratio K{a,P) of the photoionization and thermal dissociation energies as a function of 
the strength of the Coulomb interaction j3 for a large polaron bound to an impurity with different 
values of the Prohlich coupling constant a. Note that for large constants a, the calculated values 
«;(«, 0) tend to 3, in conformity with the limiting result of the adiabatic theory |^. 



FIG. 2. Visualisation of the atomic structure of the V~ center in MgO (a) and V^Mg center in 
corundum (b) as calculated by means of the INDO method. In the V~ center five nearest neighbour 
O atoms move outwards from the Mg vacancy by about 1% oq (the lattice parameter), whereas the 
O atom with a hole (shaded) is attracted to the Mg vacancy by 3% oq which makes the dominant 
contribution to the energy gain due to the lattice relaxation. 
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TABLES 

TABLE L Theoretical and experimental values of the E'th and -Eopt for a polaron bound to the 
Mg vacancy (the ^"-center) in MgO. CTAC, CTSC, INDO stand for continuum arbitrary-coupling 
theory, continuum strong-coupling theory, and microscopic quantum-chemical approach, corre- 
spondingly; experimental values are given in the last column. TiLOex is a transition energy between 
the ground and excited states. Sign "— " indicates that the respective quantities were not calcu- 
lated; indexes g and e label, correspondingly, the ground and excited states (those states are Is 
and 2p in CTAC). 



Energies (eV) 


CTAC 


INDO 


CTSC 


experiment 


EtY. 


L24 


1.6 


1.15 


1.4 i 


Eopt 


2.2 


2.2 


2.13 


2.3® 


iEopt)i 


L13 


- 


0.96 




huJex = Eopt — (-E'opt)f 


L07 


- 


1.16 




-^polariz 


0.9 


- 


0.97 




K 


L77 


1.38 


1.85 


1.64 



TABLE IL Same as Table 1 for the VMg-center in a-Al203 



Energies (eV) 


CTAC 


INDO 


experiment 


Eth 


1.22 


1.26 


0.68 |C|; 0.7 fH 


Eopt 


(2.58)* 


- 


2.58 @ 


(Kpt)f 


1.15 


- 




hiOex = Eopt — (-E'opt)f 


1.43 


- 




^polariz 


1.36 


- 




K 


2.11 


- 


3.68 



* The value of the energy Eopt 
ruh = AmQ. 



2.58ey was used to estimate the effective mass of a hole: 
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